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EFFECTS OF VIBRATIONAL NONEQUILIBRIUM 
ON THE INVISCID DESIGN OF AN AXISYMMETRIC 
NOZZLE FOR HYPERSONIC FLOW 

By John B. Anders, Jr. 

Langley Research Center 

SUMMARY 

The design of an axisymmetric, hypersonic nozzle for arc-heated air with vibrational 
nonequilibrium is computed by the vibrational sudden- freeze technique. For stagnation 
conditions near the design values, calculations based on methods for quasi- one- dimensional 
flow indicate that the chemical composition of the gas (air) freezes very early in the expan- 
sion. The vibrational mode of the molecules freezes somewhat later in the expansion. A 
method-of-characteristics calculation for equilibrium air is used to design the supersonic, 
inviscid portion of a nozzle where vibration is "suddenly frozen; " that is, upstream of the 
freeze point vibration is in equilibrium, and downstream of the freeze point vibration is 
frozen. 

This sudden-freeze nozzle is analyzed at design conditions by means of a character- 
istics calculation for frozen chemistry and vibrational nonequilibrium. The results of this 
analysis indicate that the approximate design provides reasonably uniform exit profiles at 
a Mach number about 3 percent lower than the design value of 9. The exit flow angularity 
is everywhere less than 0.1°, and the static-pressure variation across the exit is about 
10 percent of the center-line value. Analysis of nozzle performance at slightly off-design 
stagnation conditions reveals no major deterioration of the flow quality. 

INTRODUCTION 

Wind-tunnel nozzle design has commanded considerable effort over past years when 
long and laborious method- of- characteristics (MOC) calculations were performed by hand. 
As higher Mach numbers were required, two-dimensional nozzles gave way to axisymmet- 
ric designs, and approximate methods such as that described in reference 1 were used. 
With the development of high-speed digital computers, numerous techniques using the MOC 
appeared (e.g. , refs. 2 to 8). The high stagnation temperatures associated with high Mach 
number flight required calculations to account for real-gas effects (refs. 9 to 12) and for 
nonequilibrium (refs. 13 to 16). 



In order to increase the flexibility and capability of the Langley 4-foot hypersonic 
arc tunnel described in reference 17, it was desirable to replace the existing high Mach 
number conical nozzle with a lower Mach number contoured nozzle. The present report 
develops the design for the new nozzle. The high temperature of arc- heated gas indicated 
that the nozzle design must be for a real gas with possible nonequilibrium. 

One major difficulty facing the designer of a hypersonic wind-tunnel nozzle with non- 
equilibrium flow is that most available nonequilibrium MOC calculation schemes treat 
what is called the direct problem (nozzle contour is assumed given) instead of the indirect 
problem (nozzle contour is the unknown to be determined). For this reason an approxi- 
mate design scheme was employed to generate the nozzle contour, and an existing direct 
MOC program modified to include some nonequilibrium effects was used to analyze the 
nozzle flow. 

The existing tunnel complex placed many restrictions on the nozzle design. For 
example, both the arc heater and vacuum sphere limited the total mass flow of the system, 
and the operating characteristics of the arc heater limited the maximum practical stagna- 
tion pressure. The necessary trade-offs to optimize the design and the MOC calculations 
used to compute the contour are discussed. Vibrational nonequilibrium was treated in an 
approximate way, and the resulting nozzle design was analyzed by means of a direct 
vibrational- nonequilibrium MOC calculation. A boundary-layer correction to the inviscid 
contour was not computed for the nozzle, although such a correction could be quite large. 
Also, the sub sonic- nozzle contour was not considered. 

SYMBOLS 


A cross-sectional area 

a speed of sound 

Cj mass fraction of chemical species i 

Cp specific heat at constant pressure 

c p,a nondimensional specific heat, Cp ja / C p,oo 

d diameter 

e c energy per unit mass in formation of chemical species, j 

i 


2 



e f 


H 

h 

h 


4,i 

i 

M 

m 

n 


P 

P 

q 

R 

R s 

Roc 

r 

S 


s 


energy per unit mass frozen in vibration and/or chemical species 

energy per unit mass in vibration 

nondimensional vibrational energy, e^/Cp^T^ 

total enthalpy per unit mass 

static enthalpy per unit mass 

nondimensional static enthalpy, h/cp ^T^ 

enthalpy of formation for chemical species i 

index of chemical species 

Mach number 

mass flow rate 

coordinate normal to streamlines 
pressure 

nondimensional pressure, p/p^u^ 
dynamic pressure, pu^/2 
universal gas constant 
specific gas constant 

Reynolds number based on free-stream conditions at nozzle exit 
radius 

entropy per unit mass 
coordinate along streamlines 
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s nondimensional coordinate, Ar 

r 

T temperature 

T nondimensional temperature, ’T 2i /'T 00 

t time 

U velocity ratio, u/u 7 

u velocity 

u nondimensional velocity, u/uoo 

limiting velocity 

W molecular weight of mixture 

Wj molecular weight of chemical species i 

■fri chemical species production term 

nondimensional chemical species production term, w^/p^u^ 
x coordinate along nozzle axis 

x nondimensional axial coordinate, x/r* 

y coordinate normal to nozzle axis 

y nondimensional normal coordinate, y/r* 

y ratio of specific heats 

V characteristic variable 

9 flow angle 

l-i Mach angle 
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characteristic variable 


I 

p density 

t vibrational relaxation time 

Subscripts: 

a active degrees of freedom (translation and rotation) 

i chemical species i 

t stagnation condition 

v vibration 

1 arc- heated fluid 

2 room- temperature fluid 

3 arc-heated fluid mixed with room- temperature fluid 

°° free- stream conditions at nozzle exit 

Superscript: 

* throat of nozzle 


DESIGN CONSIDERATIONS 
Selection of Stagnation Conditions 

Because of the high stagnation temperatures associated with arc- heated flows, most 
test models must be water cooled. In general, the models are therefore bulky and require 
a rather large test stream to avoid blockage. In view of this, the exit diameter of the noz- 
zle was fixed at a minimum value of 15 cm. The exit Mach number was chosen to be 9 
since the arc heater could duplicate the total temperature for flight at that Mach number. 

For a real gas the area ratio required for a given Mach number is a function of the 
stagnation conditions. For heater operation at 100 atm (10.13 MN/m^) the size require- 
ment can be met by expansion to a free-stream Mach number of 9 for a mass flow of 
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approximately 0.27 kg/sec. Past experience with the present facility has shown that for 
stable operation the maximum arc-heater mass flow is about 0.19 kg/sec. At these con- 
ditions the heater can deliver 5811 kj/kg. The total mass flow of the system can be 
increased from 0.19 kg/sec by adding cold room air downstream of the arc. The final 
total enthalpy resulting from mixing the two flows will be somewhat less than 5811 kj/kg, 
depending on the cold-air mass flow. The interrelation of all these factors can be seen 
by the following simple development. From figure 1 an energy and mass balance can be 
written as 

m 1 H 1 + m 2 H 2 - rhgHg (1) 


and 


+ m 2 = m^ 


Combining equations (1) and (2) and solving for m 2 result in 


. . H 1 - H 3 

m 2 " m l 


The total mass flow is then 


m 0 = m 


H l- H 2 

1 «3 - H 2 


Now, from continuity, 


*1 H 1 - h 2 

P 3 V ' P 3 V H 3 - "2 


* m 3 
A* = d 


For axisymmetric flow, 


d*= / 4A * 


/1.275m^ Hj - H 2 
p 3 u 3 h 2 


( 2 ) 


(3) 


(4) 
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Let A/A* be the exit area ratio necessary for the desired exit Mach number at a given 
enthalpy H3. Then 


d = 



' 1.275rh 1 - H 2 A 


P 3 \* H3 " H 2 A * 


( 5 ) 


A plot of equation (5) is shown in figure 2 for an equilibrium gas expanding to 
M = 9. The largest nozzle diameter is obtained with no cold mass flow but at consider- 
able expense to the Reynolds number. Increasing m 2 increases but decreases 

Hg and d. On the basis of figure 2 the design point was chosen at Hg = 4450 kJ/kg. 
This matches the total enthalpy for M = 9 flight at an altitude of about 42.7 km. The 
exit diameter is about 17.8 cm for an equilibrium expansion with a total mass flow of 
0.25 kg/sec. If nonequilibrium effects occur, however, the exit diameter will be reduced. 
To summarize, the selected design stagnation conditions are 

T t = 3400 K 

Pj. = 100 atm 

H t - 4450 kJ/kg 

m = 0.25 kg/sec 

Moo -9 : 

At the specified stagnation conditions the composition of air is 

C M = 0.7220 
z 

C n = 0.1857 
u 2 

= 0.0702 
C Q = 0.0088 
C Ar = 0.0129 

The amount of any ionized species is negligible. 
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Determination of Gas Model 


The computer program of Lordi, Mates, and Moselle (ref. 18) for quasi-one- 
dimensional, nonequilibrium expansion provides a convenient tool for determining the 
proper gas model for the nozzle- design calculations. The program allows for 11 chemi- 
cal reactions with species N 2 , N, 0 2 , O, NO, NO + , Ar, and e". Vibration may be either 
frozen at the stagnation conditions or in equilibrium with the translational and rotational 
modes. Figure 3 shows the mass-fraction variation of four reacting species along a coni- 
cal nozzle. Vibration is assumed to be in equilibrium. The mass fractions of the three 
most abundant species N 2 , 02, and NO freeze very rapidly downstream of the nozzle throat, 
and thus the composition of the gas is essentially fixed. Figure 4 shows the energy in 
chemistry for a nonequilibrium expansion. Results are shown for nozzles with two dif- 
ferent cone angles (5° and 50°) in order to illustrate the effect of nozzle geometry in the 
nonequilibrium case. The chemical reactions freeze somewhat sooner in the 50° nozzle 
than in the 5° nozzle. 

Although the chemical reactions freeze early, even in the slower expansion of the 
5° cone, the same will not necessarily be true for vibration, since the relaxation rate, in 
general, is much faster for vibration. The program of reference 18 does not consider 
vibrational nonequilibrium, and another way of assessing this effect must be found. 

An approximate method for treating vibrational nonequilibrium is the sudden-freeze 
technique of Bray (ref, 19). This useful approximation has found many applications, and 
it was adapted in reference 20 to the nonequilibrium program of reference 18 by means of 
the freezing criterion of Phinney (ref. 21). The freeze point for vibration is defined as the 
point in the expansion where 


dey ^ ey 

"dF ~ ~ 


(6) 


In this equation, e v is the equilibrium vibrational energy of the gas, and r is the vibra- 
tional relaxation time. Usually r is determined from experimental measurements. Ref- 
erence 20 used a value for r obtained from measurements of the vibrational temperature 
(Tyj of N 2 in air for an expanding flow. For the present paper, use has been made of more 
recent electron- beam measurements of T v for N 2 in air (ref. 22). The equation that 
best fits this data is 


jp 


1.2 x 10" 10 



( 7 ) 


In this investigation equation (7) is applied to all three vibrating species, even though 0 2 , 
NO, and N 2 do not relax at the same rate. Since the vibrational-temperature measure- 
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roents are probably accurate to within ±15 percent, at best, and since does form the 
bulk of the gas (air), the use of a single vibrational relaxation rate is probably justified. 
The rates determined in shock-tube studies differ greatly from those determined in nozzle 
flows (refs. 23 and 24), and it is important that the rate used in this investigation was deter 
mined from nozzle measurements. The use of shock- tube values of r in the program of 
reference 18 results in vibrational freezing much earlier in the expansion, even upstream 
of the throat. Figure 5 shows the variation of several Tp equations, including the theo- 
retical curve of reference 25 and the experimentally determined curve of reference 26. 

No attempt has been made here to evaluate the many expressions for r. Equation (7) 
which fits the recent data of Sebacher and Guy (ref. 22), is assumed to be adequate. 

The freeze criterion of reference 21 was added to the expansion program of refer- 
ence 18, and the resulting vibrational- energy variation along a conical nozzle is shown in 
figure 6. The area ratio at which freezing occurs is obviously dependent on nozzle geom- 
etry. The chemical reactions in the expansions of figure 6 are assumed to be out of equi- 
librium, but in fact, the major species freeze quite early in the expansion as was shown in 
figure 3. Since the energy frozen in chemistry is slightly higher for the 50° expansion 
(fig. 4), a lower translational temperature occurs downstream; consequently, a lower equi- 
librium vibration curve is shown in figure 6 for the 50° nozzle. 


METHODS AND RESULTS 

f 

Method of Characteristics for Equilibrium Flow 
.• Calculating the inviscid contour of a hypersonic nozzle by the method of characteris- 

I 

tics (MOC) is a well-known technique that has been in use for many years. With the use of 
high-speed computers the procedure for equilibrium flow is rapid and straightforward. 
Nonequilibrium solutions are more difficult to obtain in practice but may be solved by the 
MOC. The MOC for an equilibrium real gas is discussed in this section, and then a more 
detailed analysis of the nonequilibrium solution is presented in a subsequent section. 

Equilibrium equations .- The equations of motion in natural coordinates for axisym- 
metric, irrotational, inviscid flow (ref. 27) are 


- 1 8U _ 96 _ sin 8 _ q 
U 9s an r 

) 

1 M _ M. = o 

U 8n 8s 


( 8 ) 
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From equation (8) the compatibility relations are 


t tan M d0- sin g . sln jjJm ■ /^ = 
U p y cos (0 ± p) 


0 


( 9 ) 


The characteristic directions are given by 

g * tan(0 ± ju) (10) 

The characteristic lines are the Mach lines. Figure 7 shows the relationship between the 
characteristic lines and the streamlines. 

Nozzle design by means of the MOC involves writing equations (9) and (10) in differ- 
ence form and specifying the necessary boundary conditions. Usually the nozzle is divided 
into several regions, as indicated in figure 8. Region IV is a region of parallel uniform 
flow at the desired Mach number, region n is a radial-flow region, and regions I and in 
are transition regions calculated by the MOC. The flow in region II is calculated by means 
of the radial- flow relation given in reference 3: 

(l-M 2 )|^+2^=0 , (11) 


This may be rewritten as 


dr 

r 


1 M 2 - 1 

2 U 


dU 


( 12 ) 


Equation (12) provides a relation between radial position and velocity and may be integrated 
directly for a perfect gas. However, for an equilibrium real gas the relation between M 
and U may not be simply expressed, and a closed-form integration cannot be performed. 

The following equation can be obtained from figure 9: 


~ = = d0\/M 2 - 1 (13) 

r tan ii v 

Combining equations (12) and (13) results in 

d (14) 
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This provides a relation between angular position and velocity. Again, for a perfect gas 
but not for an equilibrium real gas, equation (14) may be integrated immediately. Once a 
relation between M and U is supplied, the flow properties on the boundaries of the 
radial-flow region may be determined for a selected exit Mach number. In figure 8, flow 
properties along the line AB are fixed by selecting a Mach number distribution, and hence 
a U distribution. The characteristic net in regions I and III can be constructed once the 
flow properties are determined along segments AB, BC, CD, and DE. 

A computer program utilizing a MOC nozzle- calculation procedure similar to that 
just outlined was developed and reported in reference 9. The relationship between velocity 
and Mach number for an equilibrium real gas was supplied by tabulated values of l/M 2 
for various values of 1/U 2 . The Mach number distribution along line AB was assumed to 
be linear -with' slope equal to that at the beginning of the radial- flow region. The effects of 
a real gas on nozzle contours were indicated in reference 9, but the effects of frozen chem- 
istry or frozen vibration in the stagnation chamber are also of interest. Such results 
should give an idea of the importance of vibrational or chemical nonequilibrium. 

The necessary tables of l/u 2 as a function of l/M 2 for equilibrium or frozen, 
quasi- one- dimensional flows may be easily computed by means of such computer programs 
as in reference 18 or reference 28. The program of reference 18 was used to construct 
four tables assuming (1) chemical equilibrium and vibrational equilibrium, (2) chemical 
equilibrium and frozen vibration, (3) frozen chemistry and vibrational equilibrium, and 
(4) frozen chemistry and frozen vibration. Case (2) is unrealistic, since chemistry usually 
freezes before vibration, but it will serve for comparative purposes. It should be clear 
that the limiting velocity for frozen flow is reduced by the energy frozen internally, that is, 
vibrational energy or energy of reaction. The actual limiting velocity is given by 

» i - 2 ( H t - e f) 

The relation between velocity and Mach number must, therefore, vary according to the 
amount of energy frozen. This shifting effect will assume more importance in the sudden- 
freeze calculations discussed in the next section. 

Figure 10 shows the four nozzle contours calculated by means of the program of ref- 
erence 9 and the four tables described in the previous paragraph. It is immediately obvi- 
ous that freezing the vibrational energy lias a large effect on the final area ratio needed 
for Mach 9 flow. The effect seen here is more than simply the removal of a certain amount 
of energy from that available for the expansion. For vibrational equilibrium the amount of 
energy in vibration strongly influences the value of y of the gas and thus appears implic- 
itly in the relation between Mach number and area ratio. Freezing vibration in the 
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stagnation chamber for the temperature and pressure in figure 10 reduces the exit area 
ratio by approximately 50 percent for Mach 9 flow. Frozen chemistry reduces the exit 
area ratio by approximately 20 percent. 

Sudden-freeze model .- From the calculations for quasi-one- dimensional flow, it was 
apparent that the composition of the flow was essentially fixed at its stagnation- chamber 
value. This left vibration as the only nonequilibrium process. Since the sudden-freeze 
approximation of a nonequilibrium process had been successfully used in analyses for 
quasi-one- dimensional flow (refs. 20 and 29), it was adopted in the present MOC calcula- 
tions. The resulting simplifications are obvious when the compatibility relations for a 
vibrating diatomic gas are written as 



where the subscript a refers to the active degrees of freedom (translation and rotation) 
and the subscript v refers to vibration. For equilibrium or frozen flow the terms on the 
right-hand side of equations (15) are zero, and these equations reduce to equations (9). 

The sudden- freeze model, then, allows the use of the most simple form of the compatibility 
relations. The characteristic directions for this gas model remain the same as in equa- 
tions (10). The equilibrium MOC program of reference 9 is easily adapted to this tech- 
nique by simply calculating the tables of velocity as a function of Mach number for a sud- 
denly frozen gas. 

For frozen or equilibrium flow the tables of velocity as a function of Mach number 
used in the method of reference 9 are independent of geometry, but for the sudden-freeze 
model these tables are dependent on the freeze point, which, in turn, depends on the geom- 
etry. Since the actual geometry is the unknown being sought, an iterative process is nec- 
essary. As the first approximation a hyperbolic nozzle contour with 5° asymptote was 
assumed. The tables thus constructed were used in the MOC calculation to obtain a con- 
tour. This proved sufficient for the present case. However, a second iteration may be 
performed, if desired, by using the calculated contour to generate new tables and, hence, a 
new contour. 

Speed of sound .- It is useful at this point to consider the proper definition of the speed 
of sound for a nonequilibrium flow. In reference 30, Vincenti and Kruger point out that the 
speed of sound for equilibrium flow is defined only for an equilibrium state, whereas the 
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speed of sound for frozen flow is defined for an arbitrary state. For nonequilibrium flow 
there is no unique speed of propagation that depends only on the conditions in the undis- 
turbed gas. At the frozen and equilibrium limits the small- disturbance equation exhibits 
a singular drop in order, but arbitrarily close to equilibrium the frozen speed of sound 
retains its meaning, whereas the equilibrium speed of sound does not, Sedney et al. 

(ref, 31) and references 32, 33, and 34 have concluded that the frozen speed of sound is 
the correct value on which to base the characteristics, even very close to equilibrium. 

In the present case the sudden- freeze concept results in a discontinuous change in 
Mach number from the equilibrium value to the frozen value at the freeze point. In order 
to avoid this discontinuity, the frozen speed of sound was used even in the equilibrium 
region for the present MOC calculations. This seems to be consistent with the discussion 
in references 30 and 31. 

Tables of velocity as a function of Mach number .- The necessary table for use in the 
program of reference 9 may be constructed for the sudden-freeze model by using the fol- 
lowing development. The energy equation is 

H t = h + (16) 

The limiting velocity is defined by 

U; = \l2K t 


and 


U = — 


Equation (16) becomes 

1 _ H t 

u 2 " H r h 


(17) 


Now, 


h = Cp >a Ta + e v 


and 
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M 2 u£ u 1 2 

^ ~ a 2-y a R S Ta 


Therefore, 


h = 


c p,a 


_u? L_ 

M a 2 yaRs 


+ e v 


(18) 


After some algebraic manipulations, equation (17), with the use of equation (18), becomes 



( 19 ) 


which expresses the desired relationship between frozen Mach number and velocity. For 
a nonvibrating gas, equation (19) reduces to the familiar perfect-gas relation. The term 
on the right-hand side of equation (19) passes smoothly through the freeze point, since it 
contains the frozen Mach number M a . The 1/U 2 term, however, poses some difficulty 
because the limiting velocity is defined in terms of the total energy available for the expan- 
sion. At the freeze point, a certain amount of energy is ’locked" in vibration and becomes 
unavailable for conversion to velocity. Thus the limiting velocity before freezing is 

u z = ^(2Ht (20) 

and after freezing is 


u^ = pint - e f ) (21) 

Of course, if chemistry is frozen in the stagnation chamber, as has been assumed here, 
then the total stagnation enthalpy is reduced by that amount initially. 

2 i 2 

Figure 11(a) shows the variation of 1/U with 1/M a for a hyperbolic expansion 
with chemistry frozen in the stagnation chamber and vibration frozen downstream of the 
throat. Curves are shown for the two different expressions for u^ (eqs. (20) and (21)). 
Figure 11(b) shows the freeze-point region to a larger scale. The smooth curve through 
the discontinuity is strictly an arbitrarily faired transition curve. The possibilities for 
constructing the tables of velocity as a function of Mach number are 

(1) Use equation (20) for the complete expansion 

(2) Use equation (21) for the complete expansion 
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(3) Use an arbitrary transition between equations (20) and (21) in the freeze-point 
region. 

Without some rational basis for choosing the transition curve the use of method (3) seems 
unrealistic, since the choice of the transition curve amounts to assuming a vibration relax- 
ation rate. The easiest choice would be to use the equilibrium definition (eq. (20)) for the 
complete expansion, since the freeze point occurs at small values of e v . Actually, noz- 
zle designs computed by methods (1), (2), and (3) in the construction of these tables showed 
negligible differences in contour, as can be seen in figure 12. This is not surprising in 
view of the small amount of energy frozen in vibration. For example, at l/u^ = 4.0 
(fig. 11(a)) the change in M a between the two curves is only about 3 percent. At higher 
Mach numbers the percentage change is even less. 

Figure 13 shows the complete contour generated from the program of reference 9 
for a sudden-freeze gas model. As expected, the exit area ratio for this nozzle falls 
between that for the vibrational- equilibrium nozzle and that for the frozen- vibration noz- 
zle, and the exit diameter of this nozzle meets the size requirements mentioned previously. 
The design conditions fix the throat radius at 0.350 cm, and therefore the exit diameter is 
15.20 cm. Thus, vibrational freezing reduced the exit diameter by more than 14 percent 
from the original value for an equilibrium expansion (17.8 cm). Although the nozzle con- 
tour is smooth and the exit diameter is of the right order, there is no assurance that this 
approximate design will provide a Mach 9 uniform flow when vibrational- nonequilibrium 
effects are present. That is, the nonequilibrium effects on nozzle contour are only approx- 
imated by the sudden- freeze model, and consequently the contour may be wrong. In order 
to assess the usefulness of the sudden-freeze design, a detailed analysis of the effects of 
vibrational nonequilibrium on flow through a given nozzle was undertaken. The results of 
this analysis are given in the next section. 

Method of Characteristics for Nonequilibrium Flow 

The MOC for nonequilibrium nozzle flows has been treated by a number of authors 
(refs. 13 to 15 and 35 to 37, e.g.). Sedney (ref. 15) surveyed many of these methods and 
concluded that although much effort has been expended, there are still many uncertainties 
as to the particular grid scheme and starting procedures that should be used. Sedney also 
points out that most methods deal with what is called the direct problem, that is, deter- 
mining the flow in a nozzle of given contour. This is unsuitable for calculation of the noz- 
zle contour necessary to produce a specific final test- section flow (the indirect problem). 
However, the sudden-freeze nozzle design generated in the previous section can be ana- 
lyzed by means of one of these direct methods. One such direct analysis is the method of 
reference 38. This method was written for internal flows with mixing and combustion, but 
it was modified for the present problem to treat an inviscid flow with vibrational nonequi- 
librium. In order to show clearly the vibrational terms in the compatibility relations for 
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a reacting, inviscid, vibrating gas, these relations are shown in a form different from 
equations (15) but in a form similar to that presented in reference 38. 

Equations of motion .- The equations of motion are 

Energy 


p.aPU ^ + Pu2 H - -pu > Ci - ^M-i 


( 22 ) 


where 


c p,a 



i 


The additional term due to vibration is immediately evident. The other equations are 
Continuity 


1 9p 1 3u 86 sin 0 n 

p 9s + u 3s Bn y 


(23) 


n- momentum 


pu 


2 10 = _ Bp 

9S Bn 


(24) 


s- momentum 


pu 


BU _ Bp 


BS 


9S 


(25) 


State 


P = 



(26) 


Relaxation rate 


^v.i _ e v,j( T a) ~ e vq(Ty) 
9s ru 


(2V) 
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Compatibility relations .- The compatibility relations may be derived from equa- 
tions (22) to (27). Reference 39 gives a detailed derivation resulting in the following non- 
dimensional equations: 


sin p cos p a 

r dp ± dB + 

y a P 


sin 0 

y 


^-y Ci 

c T Lj v 
c p,a i a i 



. ! „ - 1 W lw: 


, a l a l 'a ^ _ q 

y (y - 1) 2 — y (y - 1) .. 2 — cos(0 ± Pa) 

'a r °° Moo pu a 00 Mf*, pu 


Equations (28) are similar to the compatibility relations given by Dash in reference 38 
except for the appearance of the vibrational- energy term. 

Design analysis .- Equations (28) were used to modify the program of reference 38 
for the case of chemically frozen air with vibrational nonequilibrium. Reference 39 
describes the changes made and gives a complete listing of the modified program. Essen- 
tially, the chemistry package in reference 38 was replaced by the vibrational terms for 
N 2 , 0 2 , and NO, so that, unlike the original program, no chemistry option was available. 

All species containing hydrogen and all diffusive terms were deleted. 

The resulting program was used to analyze the flow in the sudden-freeze nozzle 
design of figure 13. The computed exit profiles are shown in figure 14. The largest non- 
uniformity across the exit appeared in the pressure, which increased by about 10 percent 
from the center line to the inviscid-core edge. The temperature increased approximately 
3 percent away from the center line, while the velocity varied less than 0.2 percent. The 
Mach number was relatively uniform across the nozzle exit but was approximately 3 per- 
cent below the design value of 9. The dynamic-pressure variation was about 7 percent, 
and the flow angularity was negligible (less than 0.1°). For many applications the nonuni- 
formities present in this nozzle are acceptable, especially considering the approximate 
nature of its design. The vibrational- energy distribution along the nozzle center line is 
shown in figure 15. The nearly constant value of energy above p. = 40 indicates frozen 
vibration. 

Off-design analysis .- In arc- tunnel testing it is important that small variations in 
stagnation conditions do not significantly disturb the quality of the test-section flow, since, 
in many instances, exact duplication of stagnation conditions from run to run is impossible. 
A limited analysis of off- design stagnation conditions for the sudden- freeze nozzle of fig- 
ure 13 was conducted. Figure 16(a) shows the effect on the exit profiles of a 10-percent 
decrease in stagnation pressure, and figure 16(b) shows the effect of a 10-percent increase 
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in stagnation temperature. In figure 16(c) both pressure and temperature have been 
reduced by 10 percent. In all three cases flow uniformity is similar to the design case. 
Figure 16(d) shows exit profiles for a 3 5-percent- lower stagnation pressure and a 
25- percent- lower stagnation temperature. The maximum flow angularity has increased 
almost 20 percent over its maximum at the design conditions. Thus, large deviations 
from the design conditions may be harmful to the flow quality, since actual nonequilibrium 
processes are strongly dependent on stagnation conditions. 

CONCLUDING REMARKS 

The present investigation has treated the approximate design of a nozzle with non- 
equilibrium flow. Calculations for quasi- one- dimensional flow indicated that for conical 
expansions of 5° and 50° at the desired stagnation conditions of 100 atm and 3400 K, chemi- 
cal reactions were essentially frozen in the stagnation chamber and vibration was the only 
nonequilibrium process. The nonequilibrium vibration was approximated by the sudden- 
freeze concept in which the three vibrating species, N£» O 2 , an ^ NO, were all assumed to 
freeze simultaneously. An empirically determined relaxation rate for N 2 in expanding air 
was used to determine the freeze point. An existing equilibrium characteristics program 
was modified to produce a nozzle contour with vibrational sudden freezing. Since freezing 
occurred for the present conditions at low values of vibrational energy, this contour was 
found to be insensitive to shifts in the relationship between velocity and Mach number 
caused by freezing the vibrational mode. This simplified the calculations through the 
freeze-point region by allowing use of a single expression for the limiting velocity. 

The approximate contour thus generated was analyzed for design conditions by modi- 
fying an existing direct characteristics program to account for vibrational nonequilibrium. 
The exit Mach number was approximately 3 percent lower than the design value of 9 but 
was still very uniform across the nozzle exit. The temperature and velocity profiles indi- 
cated only minor variations (less than 0.2 percent in velocity and less than 3 percent in 
temperature), but the pressure varied by as much as 10 percent from the center line to the 
core edge. The dynamic-pressure variation was slightly less than 7 percent. Overall, 
however, the inviscid contour generated by this approximate method gave exit profiles con- 
sidered acceptable for the present application. 

Analysis of the operation of this nozzle at slightly off-design stagnation conditions 
(±10 percent in pressure and temperature) indicated no appreciable additional exit nonuni- 
formities. However, at a pressure 35 percent below and a temperature 25 percent below 
the design values, the maximum flow angularity was increased almost 20 percent over the 
design value. This, of course, indicates that at stagnation conditions far removed from 
the design values, the nonequilibrium effects are no longer properly accounted for in the 
contour and may cause a nonuniform exit flow. 
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This study indicates that an approximate treatment of, a single nonequilibrium proc- 
ess (vibration) can be used to compute the inviscid contour of a hypersonic, axisymmetric 
nozzle with vibrational nonequilibrium and that the contour thus generated exhibits reason- 
ably uniform exit profiles. The advantage of this approximate method lies in the fact that 
the designer may use existing computer programs and therefore save considerable time 
and effort. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., November 25, 1974. 
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Room-air 
mass flow 


Figure 1.- Schematic of energy and mass-flow balance. 




Figure 3.- Mass fractions of reacting species with vibrational equilibrium 
in a 5° conical nozzle, p^- = 100 atm; Tj- = 3400 K. 
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Figure 7.- Relation between characteristic directions and nozzle coordinate system. 



Figure 8.'- Nozzle divisions for MOC calculation. 






(b) Chemical equilibrium, frozen vibration. 

Figure 10.- Nozzle contours for equilibrium and frozen flows. 
p t = 100 atm; T t = 3400 K; = 9. 
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(b) p t = 100 atm; T t = 3740 K. 

Figure 16.- Exit profiles for off-design conditions, r* = 0.350 cm. 








